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We study a covariant formalism for the Sunyaev-Zeldovich effects developed in the previous papers 
by the present authors, and derive analytic expressions for the redistribution functions in the Thom- 
son approximation. We also explore another covariant formalism recently developed by Poutanen 
and Vurm. We show that the two formalisms are mathematically equivalent in the Thomson ap- 
proximation which is fully valid for the cosmic microwave background photon energies. The present 
finding will establish a theoretical foundation for the analysis of the Sunyaev-Zeldovich effects for 
the clusters of galaxies. 
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I. INTRODUCTION 

The Sunyaev-Zeldovich (SZ) effect (J, which arises 
from the Compton scattering of the cosmic microwave 
background (CMB) photons by hot electrons in clusters 
of galaxies (CG), provides a useful method for stud- 
ies of cosmology. For the reviews, for example, see 
BirkinshawQ and Carlstrom, Holder, and ReeseQ. The 
original SZ formula has been derived from the Kom- 
paneets equationjj] in the nonrelativistic approxima- 
tion. However, recent x-ray observations (for example, 
Schmidt et al.[5] and Allen et al.|6ij) have revealed the 
existence of high-temperature CG such as ksT e ~20keV. 
For such high-temperature CG, the relativistic correc- 
tions to the SZ effect will become extremely important. 

On the other hand, it has been known theoretically 
for some time that the relativistic corrections become 
significant at the short wave length region A < 1mm. In 
particular, the recent report on the first detection of 
the SZ effect at A < 650/im by Herschel Survey seems to 
confirm the relativistic corrections @. Therefore, reliable 
theoretical studies on the relativistic SZ effect at short 
wave length region will become extremely important for 
both existing and future observation projects. 

The relativistic SZ effect has been studied theoret- 
ically in several different approaches. Wright @ and 
Rephaeli[l(| calculated the photon frequency redistribu- 
tion function in the electron rest frame using the scat- 
tering probability derived by Chandrasekharfll], which 
is called as the radiative transfer method. The sec- 
ond approach is the relativistic generalization of the 
Kompaneets equation^, where the relativistically co- 
variant Boltzmann collisional equation is solved for the 
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photon distribution function. This approach was used 
by Challinor and LasenbyJUJ and Itoh, Kohyama, and 
Nozawa[l3j. which is called the covariant formalism. Al- 
though the two are very different approaches, the ob- 
tained results for the SZ effect agreed extremely well. 
This has been a longstanding puzzle in the field of the rel- 
ativistic study of the SZ effect for the last ten years. Very 
recently, however, Nozawa and Kohyama [l^ showed that 
the two formalisms were indeed mathematically equiva- 
lent in the Thomson approximation. This explained the 
reason why the two different approaches produced same 
results for the SZ effect even in the relativistic energies 
for electrons. 

On the other hand, there is yet another covariant 
formalism which also starts with the relativistic kinetic 
equation for photons [H|,[l6j]. Although the starting equa- 
tions are same for the two covariant formalisms, the final 
expressions for the SZ effect differ significantly from each 
other. In the present paper, we explore the two covariant 
formalisms for the thermal SZ effect in detail. We will 
show that the two formalisms are indeed mathematically 
equivalent in the Thomson limit, which is fully valid for 
the CMB photon energies. It is also found that the two 
formalisms give the same expression for the kinematical 
SZ effect in the Thomson limit. This will conclude that 
the existing formalisms 0, [3 [H| for the SZ effects for the 
CG are equivalent in the Thomson approximation. The 
present finding will establish a theoretical foundation for 
the analysis of the SZ effects for the CG. 

The present paper is organized as follows: In Sec. II, we 
derive analytic expressions of the redistribution functions 
for the SZ effects with the covariant formalism of Nozawa, 
Kohyama, and Itoh[l3, [13 ■ In Sec. Ill, we then rewrite 
the redistribution functions of Poutanen and Vurmflrij 
formalism in the Thomson approximation. We will show 
that the redistribution functions in the two formalisms 
are identical. The isotropic photon scattering approx- 
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imation is also presented. Finally, concluding remarks 
are given in Sec. IV. 



II. COVARIANT FORMALISM OF NOZAWA, 
KOHYAMA, AND ITOH 

A. Rate equations in the Thomson approximation 

Recently, the present authors have published a series 
of papers [2 QjlllSlliil on ^ ne Compton scattering of the 
CMB photons. In the present subsection, we summarize 
the formalism (denoted NKI formalism hereafter) for the 
SZ effects of the CG. 

In Nozawa and Kohvamaflij. it was shown that the 
covariant formalism [l3| and radiative transfer methodQ 
were mathematically equivalent in the following (Thom- 
son) approximation: 



UJ 

7— < 1 

m 
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7 = 



(1) 
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where w is the photon energy, 7 is the Lorentz factor, and 
(3 and m are the velocity and rest mass of the electron, 
respectively. Throughout this paper, we use the natural 
unit H = c = 1, unless otherwise stated explicitly. For 
the CMB photons, Eq. ([I) is fully valid from nonrela- 
tivistic electrons to extreme-relativistic electrons of the 
order of TeV region. The kinematics in the present pa- 
per are defined as follows: The CG is moving with a bulk 
velocity (3c (=vc/c) with respect to the CMB frame. As 
a reference system, we choose the system that is fixed 
to the CMB. The z axis is fixed to a line connecting the 
observer and the center of mass of the CG. (We assume 
that the observer is fixed to the CMB frame.) We choose 
the positive direction of the z axis as the direction from 
the observer to the CG. 

The rate equation for the photon distribution function 
n(x) was derived by the present authors [3, H3] under 
the assumption of Eq. (TTJ). Here, x — w/fcsTcMB is the 
photon energy in units of the thermal energy of the CMB, 
and s is the frequency shift defined by e s = x'/x. We 
recall the results here to make the present paper more 
self-contained. They are given as follows UJ: 



dn(x) 



dsP 1 {s,p c ,z)[n{e s x)-n{x)} , (3) 



dr = n e aTdt , 



(4) 
(5) 



where f3c,z is the bulk velocity of the CG parallel to 
the z-axis, n e is the electron number density, and ot 
is the Thomson scattering cross section. It should be 
noted that 0(j3 c ) and higher-order contributions were 
neglected in deriving Eq. (0} , because f3c "C 1 is satisfied 
for most of the CG. 



In Eq. (UJ, Pl(s) is the frequency redistribution func- 
tion, and Pi k( s ) is the term which appears in the case 
of non-zero bulk motions. They are defined by 



P 1 (s)= / d/?/3 2 7 5 pe(7)P(s,/3), 



Pi,k(s) = / 

JB-rr, 



d/?/3 2 7 V(7)PK(s,/3) : 



(6) 
(7) 



where 



2/V ' 



1 



Pk(s,/3) = <K/3) \p K (8,p)-P(8,P) 



PkM) 



2/3 7 6 



/(Mo.Mo) = 



8 



Ms) 1 

ft ( s ) ^°(l-/3^o) 3 
1 



i + moMo' + ^-^oXi-mS 2 ) 



r/(MoX), (8) 
(9) 

/(/xo, Mo) , (10) 
(11) 



In Eq. ©, 6(13) is a factor related to the electron dis- 
tribution function, which is, in general, a function of 0. 
The explicit forms are given by Nozawa and Kohvama[14j 
for three different electron distribution functions. In the 
present paper, for simplicity, we use the thermal electron 
distribution for p e (j), which implies 



6({3) 



7 



(12) 



where 6 e = kBT e /mc 2 is the electron thermal energy of 
the CG in units of the electron rest energy. The electron 
thermal distribution function is 



Me (7) 



1 



,K 2 (l/9 e ) 



exp(-j/0 e ) , 



(13) 



where K 2 (z) is the modified Bessel function of the second 
kind. Variables appearing in Eqs. (|6|) - (|11[) are summa- 
rized as follows: 



/3 mm = (l-e Hs| )/(l + e- |s| ): 

M = [l-e s (l-/3Mo)]//?, 
-1 

[l-e- s (l + (3))/fi 



«hh|'" 6 "' (1 " w 



for s < 
for s > 

for s < 
for s > 



(14) 
(15) 

(16) 
(17) 



For most of CG , t -C 1 is satisfied. Then one obtains 
the following solution for Eq. (j3j): 



An(x) = n(x) — n (x) 

= An t (x) + /3 c , z An k (x) , 



An t (x) = t 
Ank(x) = t 



dsPi(s)n (e s x) — n (x) 
dsPi.K(s)no(e s x) — uq(x 



(18) 

(19) 
(20) 
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where 



n (x) 



e x - 1 



(21) 



is the initial Planckian photon distribution function. 
Equations (fT9f and (|20|) correspond to the thermal SZ 
effect and kinematical SZ effect, respectively. 



B. Analytic expressions for P(s,f3) and Pk{s,P) 

First, we show the analytic expression for the redistri- 
bution function P(s, ft) in the NKI formalism. In Eq. ©, 
the integral of fXo can be done analytically. The explicit 
form was derived by Enfilin and Biermann[l8| and also 
given by Eqs. (18) and (19) in Ref . [l9| . In the present 
paper, we rewrite the expression into the following form: 



P(s,ft) = e 3s / 2 Q(s,ft), 
3 



(22) 



Q(s,p) 



16/3 2 7 4 
4 



m< 3 + rP 



2/3 



1-/3 



-ij|-lsinh^+2(/3 2 -3)| S |coshi 



+ (9 - /3 2 + 4 7 2 /3 4 ) sinh 



(23) 



Note that Q(s, /3) is an even function on s, namely 
Q(s,ft) = Q(—s,ft). Thus, the non-symmetric structure 
of the redistribution function Pi (s) is due to the function 
e 3s/2 m gq_ (J^|) , because /3 mi „ in Eq. © is also an even 
function on s. 

Similarly, we show the analytic expression for the re- 
distribution function Pk(s, /3) in the NKI formalism. The 
explicit form for Pk(s, /3) was given by Eqs. (25) and (26) 
in Ref. [13]. In the present paper, we rewrite the expres- 
sion into the following form: 



PkiaJ) = e 2s Q K (s,ft), 



(24) 



16/3 2 7 4 



2 

w 

1 



1-/3 
2, -ft 2 



2/3 7 2 



^ i -3 + 3/3 2 + 2 7 2 ^ 4 + ^r^r- In 



1 + 



2/3 7 2 



41 



1-/3J 



cosh ; 



(I 1 7 2 1 



3-/3 2 
7 2 

32 



I cosh i 



(9 - 10/3 2 - /3 4 - 4 7 2 /3 6 ) sinh |s| 



(25) 



Note that Qk( s j/3) is also an even function on s. In- 
serting Eqs. dll|) and (JUJ) into Eq. ©, one obtains the 
analylic expression for Pk(s,/3). Thus, the full analytic 



expressions for P(s, /3) and Pk(s, ft) have been derived in 
the NKI formalism. 

Before closing the present section, it is worth to men- 
tion the following. The thermal SZ effect (Eq. (fT9"]l) and 
kinematical SZ effect (Eq. ([20])) are obtained by the dou- 
ble integrals over the variables s and /3, because the an- 
alytic forms for P(s, ft) and P^s^ft) are given in terms 
of Eqs. (|22[) and (|24p . This makes the numerical calcula- 
tions of the SZ effects extremely fast, which will be quite 
useful for the analysis of the observation data. The nu- 
merical programs are available upon request from one of 
the present authors (S. N.). 



III. EQUVALENCE BETWEEN TWO 
FORMALISMS IN THE THOMSON 
APPROXIMATION 

A. Thermal SZ effect in the Thomson 
approximation 

In the present subsection, we show that the formalism 
shown by Poutanen and Vurm|16| (denoted PV hereafter) 
for the thermal SZ effect (ftc = case) is equivalent in 
the Thomson approximation to the NKI formalism. Be- 
fore to proceed the calculation, it should be remarked the 
following. In the present paper, the photon energies are 
expressed in units of the thermal energy of the CMB. On 
the other hand, all energy variables in the PV paper are 
in units of the electron rest energy. In order to make the 
present paper self-consistent, we introduce new variables 
xpv and x'py which correspond to x and x\ in the PV 
paper as follows: 



x PV = 



mc* 



XOCMB 



= x'OcMB 



(26) 
(27) 



where 6» CM b = k B T CMB /mc 2 . 

The source function is defined by Eq. (174) of the PV 
paper. We rewrite the expression with variables defined 
in the present paper as follows: 



S P y(x) 



Ppv(s,ft) = 



ds Pi^y(s)n (e s x) . 



dftft^ b Pe (j)P PV (s,ft), 



i6ftr 



:Xpye 



(28) 



(29) 



(30) 



where /i m i n and /i max are the minimum and maximum 
values of the cosine of the scattering angle, and the 
function Rq is given by Eq. (E2) of the PV paper. It 
should be noted that the angular integrals in Eq. (|30|) 
were performed by Brinkmann [2ll| and also by Nagirner 
and Poutanen(22|. The source function is related to the 
change of the photon occupation number by 



An t ,pv{x) = Spy(x) - rn Q (x) , 



(31) 
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which corresponds to Eq. ([19)) . 

In the Thomson approximation (xp V < 1), one finds 

Mmin = -1 , (32) 
e~ s (e s - l) 2 

The leading-order terms (1/xpy terms) of Ro can be 
rewritten in the Thomson approximation by 

R^l~--,+Jr-A^-?>Q 2 ) ■ (34) 



Q q a 3 2q 2 a 5 \ a 2 

Note that a± « a and d± w d were used in deriving 
Eq. (|34l) . The variables appearing in Eq. (|34p are ex- 
pressed by: 



Q = xpvVl + e 2s 
d = x P v7(e s + 1) , 
a = \/7 2 + r , 
q = xl Y e s (l - /i) , 

1 -/! 



2e s /i, 



(35) 
(36) 
(37) 
(38) 

(39) 



Inserting Eq. (|34[) into Eq. (|30p and expressing in terms 
of the variables of Eqs. (f3"5)) - (15i?l) . one can write the re- 
distribution function as follows: 

[4/3 7 2 e 2s /i - 4/3 7 3 e s (e s + 1)I 2 



32^ 2 7 4 . 

+^{e s + l) 3 / 3 - h\e s - 1)V + l)h 



2/3 7 V(e s + 1)7 5 



where 1-^ , . . . , I 5 are the integrals defined by 



h 

h 
h 
h 
h 



dfi 
dfx 



1 



(l + e 2s -2e s ^) 1 / 2 ' 
1 



, (l-,i)( 7 2 +r) 3 / 2 
1 



dfi 
dfi 
dfi 



(1 


-/i) 2 ( 7 2 


f r) 7 /2 




1 




(1 


-A*) 2 (7 2 


f r )5/2 




1 




(1 


- M)(7 2 4 


r )5/2 • 



(40) 

(41) 
(42) 
(43) 
(44) 
(45) 



Their elementary integrals can be performed. One has as 
follows: 



h 



[{e s + l)fi-\e s , 



h — hi + I22 , 

_ 2 1 e- s (e s ~ \) 2 \e s - 1| 
21 ~ 7 ~ 2/3V (e s + 1) ' 

I22 = + ^e- s (e s + l)\e s 1| 



/3 2 7 2^ 3 7 3 
1 



(46) 

(47) 
(48) 



(49) 



1 1 (e s - l) 4 |e s - 1| 
7 5 /? 5 7 5 (e s + l) 5 



1 (e s - l) 2 \e s - l\ 

33 -,3 



7 3 ^3 ( e « + l)3 
^5 = ^51 4 I52 , 

2 1 e- s (e s - l) 4 \e s - 1| 



/51 
I52 



7 3 2/3 5 7 5 (e s + 1) 3 

"^V (3 " 2/3 ^ + 05 7 5 e .s + 1 
1 

3 

/3V 



(50) 

(51) 

(52) 
(53) 



477^-( e s + l)| e s -l| 



where 



Xa = In 



(54) 



(55) 



Inserting Eqs. (|4"6"j) - (f54"|) into Eq. (|4"0")) . and performing 
a straightforward calculation, one finds that Ppv(s,/3) 
coincides with P(s, (3) given by Eqs. ([22]) and <j23j). Thus, 
the PV formalism is equivalent to the NKI formalism in 
the Thomson approximation for the thermal SZ effect. It 
should be remarked that ^3, I4 and /51 terms in Eq. (|40"]) 
have cancelled out completely in deriving Ppy(s,f3). 

Before closing the present subsection, it should re- 
marked that the equivalence of the two formalisms is 
valid not only for electrons in the thermal distribution 
but also in other distributions such as power-laws. 



B. Kinematical SZ effect in the Thomson 
approximation 

In the present subsection, we show that the PV 
formalism[16] for the kinematical SZ effect (/3c ^ case) 
is also equivalent to the NKI formalism. In the PV paper, 
0(/3c) and 0(/3 c ) terms are calculated for the kinemat- 
ical SZ effect. On the other hand, as discussed in the 
previous section, (3c <C 1 is satisfied for most of the CG. 
In Nozawa, Itoh, and Kohyama[23], it was shown that 
0(/3 c ) contribution can be safely neglected. Therefore, 
we restrict ourselves to the 0(/?c) terms. 

The source function of 0(/3 c ) in Eq. (174) of the PV 
paper is expressed by 



Sk,pv(x) = (3c, z t I dsPi, KPV (s)n (e s x) , (56) 

O 

d/3/3 2 7 5 p e (7)^K,Pv(s,/3), (57) 



-Pl,KPv(s) = / 



P K ,pv(s,/3) = 7r fk, PV (s,/3)-fpv(s,/3) 



(58) 



32^4 V e 3s - , (59) 
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where Ppy(s,(3) is given by Eq. (l30|) . and the function 
i?s is given by Eq. (E3) of the PV paper. It should be 
remarked that (3c, * = — ^6 was used in deriving Eq. (J56j) , 
where 77 is the cosine of the polar angle of the initial 
photon and (3b = (3c- The source function is related to 
the photon occupation number by 

fc^An^pvta;) = Sk,pv(x) - (3c,zTn (x) , (60) 

which corresponds to Eq. ([20)) 



The leading-order terms (l/xp v terms) of Rs can be 
rewritten in the Thomson approximation by 



R ~ M 4 



qa 



M 2 
q 2 a 5 



q 2 a 3 



(61) 



Inserting Eq. (I61[) into Eq. (|59[) and expressing in terms 
of the variables of Eqs. (|35)) -(f39 |) . one can write the re- 
distribution function as follows: 



Pk,pv(s,(3) 



3e s 



32/?V 



4/3 7 V s (e s + l)i^i - 4(3je s K 2 



+/3 7 3 (e s + 1) 2 K 3 -P 1 (e s ~1) 2 K 4 



where K±. 
K x = 

K 2 = 
K 3 = 
K 4 = 
K 5 = 



-2(3 ie 8 K 5 



. , K$ are the integrals defined by 



-1 



dfi 



dfi 



1 



(1 



2e s /x) 3 / 2 



(7 3 + r) 1 /a ' 
1 



dfi- 



-1 

Cms 



1 



(/// 



m)(7 2 
1 



„)5/2 



(l-/i)( 7 2 +r) 3 /2 ' 
1 



( 7 2 + r) 3 / 2 



(63) 
(64) 
(65) 
(66) 
(67) 



Their elementary integrals can be performed. One has 
as follows: 



Kx = 



2(3j 2 
K 2 = —I22 , 
K 3 = h, 
Ki = h , 
K 5 = -hi . 



(l- 27 2 )|e s -l|+2/3 7 2 



1 



(68) 

(69) 
(70) 
(71) 
(72) 



Inserting Eqs. (|6"5|) (|72p into Eq. (R)2"j) . and performing a 
straightforward calculation, one finds that Pk,pv(s, (3) 
coincides with P K (s,(3) given by Eqs. (JMJ) and (|2"5|) . 
Therefore 



Pk(s,(3) = j P Kj pv(s,/3). 



(73) 



Thus, the PV formalism is also equivalent to the NKI 
formalism in the Thomson approximation for the kine- 
matical SZ effect. 



C. Redistribution functions in the isotropic 
scattering approximation 

The isotropic scattering (ISO) approximation f or p ho- 
tons has been studied, for example, in Rephaeli[10| for 
the thermal SZ effect. In the NKI formalism, the ISO 
approximation can be imposed by averaging over the ini- 
tial photon solid angle, namely 



If 1 



(74) 



where /(mo,A*o) is defined by Eq. ([TT]). Thus, the ISO 
approximation in the NKI formalism is obtained by re- 
placing 



/(Mo,Mo) 



(75) 



Inserting Eq. (f75|) into Eqs. © and (ITU)) , one obtains 
the redistribution functions in the ISO approximation as 
follows: 



(62) P iso (s,(3) 



e 3s / 2 
2/?V 

„2s 



S \s\ 
(3 cosh - — sinh — 



(76) 



PK, iso (s,(3) 



^ [2,3 coshs - (1 + /3 2 ) sinh H] ,(77) 



K.iso 



Pi 



K.iso 



(s,(3)-P iso (s,(3) 



(78) 



On the other hand, the expressions of the ISO approx- 
imation in the PV formalism are given by Eqs. (145) and 
(146) of the PV paper. In the Thomson approximation, 
they can be rewritten by 



Rq 



3Q ' 
8d 
3Q 3 



(79) 
(80) 



Inserting Eqs. ([79]) and (J80J) into Eq. ([30]) and ([59j) and 
repeating the same straightforward calculations done in 
the previous subsections, one finally obtains the same 
expressions as Eqs. (f76]) and ([77]) . 

Finally, we study the accuracy of the ISO approxima- 
tion compared with the full calculation. We define the 
change of the spectral intensity function AI(x) by 



AI(x) = M t (x) + Pc,.*I k (x) , 
AI t (x)=Iox 3 An t (x), 
AI k (x) = I Q x 3 An k (x) , 



(81) 
(82) 
(83) 



where Iq = (fcsTcMB) 3 /27r 2 , and Ant(x) and An k (x) are 
given by Eqs. (fT9]) and ([20]) . respectively. 

In Fig. 1, we plot AI t (x) and AI k (x) for a typical 
electron temperature fcsT e =15keV of the CG. In Fig. la, 
the solid curve is the full calculation with Eq. (|2"21 for 
P(s,/3), and the dashed curve is the ISO approximation 
with Eq. (|T6"]) for P(s,(3). In Fig. lb, the solid curve is 
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FIG. 1: Dependences of AI t {x) and AI k (x) for k B T e = 15keV 
as a function of x. Figs, la and lb are AIt(x) and Alk(x), 
respectively. The solid curve is the full calculation, and the 
dashed curve is the ISO approximation. 

the full calculation with Eq. ([24]) for Pk(s,/3), and the 
dashed curve is the ISO approximation with Eq. (|77)) for 
Pk(s, /3). It can be seen that the ISO approximation is an 
excellent approximation for the SZ effects of the CG. Our 
finding confirms the results of the PV paper on the ISO 
approximation. The errors of the ISO approximation are 



1% at the peak positions. 



IV. CONCLUDING REMARKS 

We studied the covariant formalism for the SZ effects 
developed by the present authors (NKI). We obtained 
the analytic expressions for the redistribution functions 
P(s,/3) and Pk(s,(3) of the thermal and kinematical SZ 
effects for the CG, respectively. In the NKI formalism, 
the change of the photon occupation number An(x) and 
the spectral intensity function AI(x) can be expressed 
by the double integral forms. Thus, their numerical cal- 
culations became extremely fast compared with previous 
numerical calculations. The numerical programs will be 
quite useful for the analysis of existing and future obser- 
vation data of the SZ effects. They are available upon 
request from one of the present authors (S. N.). 

We also studied another covariant formalism for the 
SZ effects recently developed by Poutanen and Vurm. 
We showed that the two formalisms were mathematically 
identical in the Thomson approximation, which is fully 
valid for the CMB photon energies. This has concluded 
that the existing formalisms [3, Gil °f the SZ effects 
for the CG are equivalent in the Thomson approxima- 
tion. Thus, the present finding establishes a theoretical 
foundation for the analysis of the SZ effects for the CG. 

Finally, we calculated the SZ effects in the isotropic 
photon scattering (ISO) approximation. The change of 
the spectral intensity function AI(x) was calculated both 
in the full expression and in the ISO approximation. It 
was shown that the ISO approximation was an excellent 
approximation compared with the full calculation. The 
errors were 1% for AI t {x) and Alk(x) at the peak posi- 
tions. 
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